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Abstract. We present interpretations of the line-of-sight velocity distri- 
£NJ ■ bution of stars in the Draco dwarf spheroidal galaxy in terms of CDM 

^ ' and MOND assuming constant mass-to-light ratio and anisotropy. We 

estimate the two parameters by fitting the line-of-sight velocity dis- 
persion and kurtosis profiles for stellar samples differing by a number 
of stars rejected as interlopers. The results of the fitting procedure 
for CDM, high mass-to-light ratio (131-14f solar units in V-band) and 
weakly tangential orbits, are similar for different samples, but the qual- 
ity of the fit is much worse when fewer interlopers are removed. For 
MOND, the derived mass-to-light ratio (2f solar units) is too large to 
be explained by the stellar content of the galaxy. 
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£3 . 1 Introduction 

r/a ■ 

• • , The Draco dwarf is a generic example of the class of dwarf spheroidal (dSph) 
. ' S a l ax i es °f the Local Group and as such has been a subject of intensive study in 
, recent years. The object is interesting from the point of view of theories of structure 
I_h ' formation due to its large dark matter content and its implications for the missing 
satellites problem. Dwarf spheroidals are also in the regime of low accelerations and 
therefore can provide critical tests for the alternatives to cold dark matter (CDM) 
scenarios, such as the Modified Newtonian Dynamics (MOND). New observations 
of the object performed recently provided a larger kinematic sample of stellar 
velocities (Kleyna et al. 2002; Wilkinson et al. 2004) and a better determination 
of its luminosity profile and shape (Odenkirchen et al. 2001). 

This has encouraged us to study the Draco dwarf using a method of veloc- 
ity moments (Lokas, Mamon & Prada 2005) first applied to the Coma cluster of 
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Fig. 1. Left panel: line-of-sight velocities versus distances for 207 stars in the sample of 
Wilkinson et al. (2004) with lines separating supposed interlopers from members. Right 
panel: velocity moments calculated from different samples with 0, 4 and 18 interlopers 
removed (solid, dashed and dotted lines respectively). 

galaxies by Lokas & Mamon (2003) and tested with iV-body simulations by San- 
chis, Lokas & Mamon (2004) and Wojtak et al. (2005). The method relies on 
joint fitting the velocity dispersion and kurtosis profiles which allows us to break 
the degeneracy between the mass distribution and velocity anisotropy. We have 
assumed that the dark matter density profile is given by a formula with an inner 
cusp and exponential cut-off at large distances recently proposed by Kazantzidis 
et al. (2004) as a result of their simulations of tidal stripping of dwarfs by the po- 
tential of a Milky Way size galaxy. We found that the results depend dramatically 
on the sample of stars under consideration and concluded that the larger sam- 
ples are probably contaminated by unbound stars because the velocity moments 
constructed from them arc inconsistent. 

Since this dark matter profile is not very different in shape from the distribution 
of stars given by the Sersic profile, in this contribution we reanalyze the same data 
assuming a simpler hypothesis that mass traces light and using different binning. 
This case is directly comparable to the alternative interpretation of the data in 
terms of MOND (Sanders & McGaugh 2004) which also assumes constant mass- 
to-light ratio but modifies the gravitational acceleration. We examine the two 
hypotheses subsequently in the next two sections. 

2 Cold Dark Matter 

The kinematic data we have used are plotted in the left panel of Fig. ^ which 
shows the line-of-sight velocities versus distances of 207 stars counted as members 
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Table 1. Fitted parameters in the case of CDM 
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Fig. 2. Velocity moments for the most restrictive sample (with 18 interlopers removed) 
together with the best-fitting model in the case of CDM. 



of Draco by Wilkinson et al. (2004). Since some of the stars clearly are discrepant 
from the main body of the galaxy we further exclude some of them as interlopers. 
The selection is indicated by the solid and dashed lines intended to follow the 
overall shape of such diagrams for gravitationally bound object. The line-of-sight 
velocity moments, dispersion <j\ os and kurtosis k\ OSi calculated from the different 
samples thus obtained are plotted in the right panel of the Figure. In each case 
we divide the data into 6 radial bins with 30 stars. We see that the results in the 
outer bins from which the interlopers are removed are affected dramatically and 
both moments are significantly reduced in value. We have assumed that the binary 
fraction in the stellar sample is small and its effect on the moment negligible in 
comparison with other uncertainties (but see the discussion in Lokas et al. 2005). 

Assigning sampling errors to the moments we fit them with the models based on 
solutions of the Jeans equations characterized by only two constant parameters, the 
mass-to-light ratio in V-band, M/Lv, and the anisotropy parameter (3 = l—dg/a 2 .. 
The remaining assumptions and adopted parameters are as in Lokas et al. (2005). 
The results for different samples are listed in Tabled 
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Fig. 3. The la, 2a and 3cr confidence regions in the parameter plane for the most 
restrictive sample (with 18 interlopers removed) obtained from fitting velocity dispersion 
alone (left panel) and both moments (right panel) in the case of CDM. 



Fitting first the velocity dispersion alone we find a stronger preference for tan- 
gential orbits in the case of samples with smaller number of interlopers removed. 
This is understandable: since the shape of the mass profile is constrained to be 
strongly decreasing by the assumption of constant MjLy the strongly increasing 
dispersion profile induces tangential orbits. The tendency disappears, however, 
once the kurtosis is included in the analysis; then all samples yield similar best- 
fitting parameters although the quality of the fit is much worse for the less restric- 
tive samples. The reason for this is the fact that for these samples the moments 
seem inconsistent: if the increasing dispersion profile is due to tangential orbits 
then this should result in decreasing instead of increasing kurtosis profile (see the 
discussion in Lokas et al. 2005). 

The velocity moments together with the best-fitting models for the case of the 
most restrictive sample (with 18 interlopers) are shown in Fig. El Fig-Glillustrates 
the benefit from including the kurtosis in the analysis for the same sample. While 
the best-fitting parameters are similar in both cases, the confidence limits for (3 
one can read from the contours shown in the Figure are much more constrained 
from fitting both moments (right panel) than from fitting the dispersion alone (left 
panel). 

The mass-to-light ratio MjLy = 131M Q / Lq found here for the most restric- 
tive sample is similar to the minimum value of 134 reached at 10 arcmin for the 
corresponding sample of 189 stars considered by Lokas et al. (2005) (see their Fig. 
6). Note however, that the quality of the fit is somewhat worse here because of a 
smaller number of parameters. 
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Fig. 4. The la, 2<r and 3cr confidence regions in the parameter plane for the most 
restrictive sample (with 18 interlopers removed) obtained from fitting velocity dispersion 
in the case of MOND assuming standard M/Lv (left panel) and ao (right panel). 



3 Modified Newtonian Dynamics 



The lowest order Jeans equation which models the velocity dispersion can be sim- 
ply modified in the case of MOND by replacing the Newtonian gravitational ac- 
celeration (?n with <7m related by #n = Ai(sM/ao)<?M where the function /i is chosen 
as fi(x) — x(l + a; 2 ) -1 / 2 (Milgrom 1983) and a is the characteristic acceleration 
scale of MOND supposed to be universal. The value found to fit well most rota- 
tion curve data for spiral galaxies is ao = 1.2 x 10~ 8 cm s~ 2 (Sanders & McGaugh 
2004). 

We have fitted the velocity dispersion data for the most restrictive sample, first 
assuming M/Ly = 3M0/L0, a value characteristic of the Draco stellar sample, 
and finding the best-fitting parameters ao and (3. Next, we adopted the canonical 
value ao = 1.2 x 10~ 8 cm s~ 2 and fitted M/Ly and (3. The results, in terms of 
best-fitting parameters and confidence regions are shown in Fig.fJJ The quality of 
the fit, x 2 /N — 9.1/4, is the same for both cases. 

In agreement with the earlier findings of Lokas (2001, 2002) using older kine- 
matic data, we conclude that the velocity dispersion profile of Draco cannot be 
well reproduced with reasonable values of either ao or M/Ly. The best- fitting 
ao = 8.7 x 10~ 8 cm s -2 and MjLy = 21.4 M /L are about an order of magni- 
tude larger than expected and the confidence regions are rather narrow excluding 
the standard values at more than 3cr level. These values are even higher (while 
orbits more tangential and quality of the fits worse) for the other samples with 
less interlopers removed. 
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4 Conclusion 

We conclude from this study that the constant mass-to-light ratio models provide 
a reasonably good fit (with x 2 /N — 15.4/10) to the kinematic data for Draco, 
but only for the most restrictive sample (with 18 interlopers removed). For other 
samples, although similar best-fitting parameters are found from the joint fitting 
of both moments, the fits are much worse (x 2 /N = 27.6/10 and x 2 /N = 39.1/10 
respectively for samples with 4 and interlopers removed) . This result is in quali- 
tative agreement with the conclusion reached by Lokas et al. (2005) who considered 
a wider class of dark matter distributions. 

The velocity dispersion data are also poorly fitted by MOND. The best- fitting 
parameters for the most restrictive sample, either ao or M/Ly are found to be 
about an order of magnitude larger than expected. However, before we conclude 
that the case of Draco really falsifies MOND other possibilities have to be consid- 
ered. The present results might mean that ao is not really a universal constant 
in MOND as previously claimed and takes different values for different classes of 
objects. Another possibility to be explored is that in modifying the gravitational 
potential according to MOND, the influence of the Milky Way (here neglected) has 
to be taken into account. In addition, since in MOND the mass is proportional to 
the fourth power of velocity dispersion any interlopers still present in our sample 
may artificially inflate the mass more strongly than in the case of CDM. 

The most interesting issue to address in the future research on dSph galaxies 
is the determination of their dark matter content without the uncertainties due 
to unbound stars. This could be done either by modelling more distant dwarfs 
which are not under direct influence of the giant galaxies of the Local Group or 
by careful removal of interlopers aided perhaps by simulations of their origin. 
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